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1. Introduction 

With the help of the multiphase chemical evolution model (MoUa, 
Hardy & Beauchamps 1999) we have derived the evolutionary histories 
of the galaxies NGC 4303, NGC 4321 y NGC 4535. With these histories 
and an evolutionary synthesis model, we were able to reproduce their 
observed radial distributions of the spectral indices Mg2 and Fe5270. 

Chemical evolution models however may fit well the present day ob- 
servational characteristics of galaxies without discriminating among very 
different star formation histories. Furthermore, the spectral indices Mg2 
and Fe5270 exhibit the age-metallicity degeneracy problem, implying 
that discriminating ages from metallicity with only these indices in sin- 
gle stellar populations might be impossible. Our objective here is to 
make sure that the implied evolutionary histories represent well these 
galaxies. 

2. The multiphase model 

The multiphase scenario begins with a gaseous protogalaxy whose 
mass is calculated from a rotation curve. The gas collapses onto the 
equatorial plane, thus forming the disk, at a rate which depends on the 
total mass: = [MMWGn/2 

The gas collapses more rapidly in the central regions than in the outer 
disk, so we assume: TcouiJ^) = tq exp ((R — Ro)/Ad). Thus, we must 
select two parameters, r^^gai and Ad, for a particular galaxy model. 



* The National Radio Astronomy Observatory is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. 
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Stars form in the halo from the primordial gas following a Schmidt 
law. In the disk, molecular clouds form from the diffuse gas while stars 
form in a second step as a result of cloud-cloud collisions. Stars also 
form via the interaction of massive stars with molecular clouds. The 
former is a local process governed by an efficiency constant valid for 
all galaxies. The same is valid for the halo star formation efficiency. 
Thus, for a given galaxy only the two efficiencies, namely and sh, 
must be chosen in order to compute a model. They are larger for earlier 
morphological types than for the later ones, being well represented by 
a probability function depending on the Hubble Type T: = e~^^/^^ 
and en = e~^^/^ (see MoUa, Di'az k. Ferrini 2002 for details). 

Summarizing, to run a model for a galaxy, i.e NGC 4303, we need 
3 input parameters: 1) The collapse time scale tq, determined by the 
total mass obtained from V(R); 2) The Hubble type T, with which we 
determine the set of efficiencies (e^, en)', 3) The scale length A^- 

3. The uniqueness problem of models for NGC 
4303 

The problem is that uncertainties exist in the selection of the above 
parameters. For the galaxy NGC 4303, the optical and radio rotation 
curves look different: the maximum rotation velocity is a factor 1.4 larger 
in the second case. This in turn implies that the collapse time scale may 
be different than the one selected 6 Gyr) in the first case. It follows 
that the efficiencies will also have to change in order to fit the present 
day observational data. Thus, alternative models would seem possible. 

We have varied these input parameters within a reasonable range cen- 
tered around those values selected in Molla et al. (1999) for NGC 4303. 
Models with all possible combinations of tq =1, 4, 8, 12 and 16 Gyr, 
A/p=l,4,8,12,and 16 kpc, and T from 1 to 10 with a variation step of 0.5, 
that is 20 possible values for the set (e^, en), have been ran, thus doing 
a total of 500 models. (See details in MoUa & Hardy 2002). 

4. The method 

In order to select the best models out of the above 500, we compare 
the radial distributions for the oxygen, diffuse and molecular gas, and 
star formation rate resulting form our models with the observations. 

This comparison is performed using the statistical indicator x^, which 
measures the goodness of fit via the differences between the model results 
and the data, taking into account the measured dispersion of data. We 
show in Table 1 the characteristics and parameters of models falling 
within the 97.5 % confidence region (i.e., a 90% of the combined values). 



Determining the star formation histories in disks 3 



Table 1. Parameters of the Selected Models (P > 90%) 
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We display in column (1) the number identifier for the model. The 
characteristics collapse time scale, tq, and the scale length, Xd, are in 
columns (2), and (3) respectively. Column (4) contains the value of T. 
The and en efficiencies are in columns (5) and (6). 

The probabilities of these distributions to be in a region around the 
minimum value of x^, for each one of our observational constraints, 
are listed in columns (7) to (10). By selecting the best models as those 
falling within the 90% confidence region we reduce the number of possible 
models to 19. All other models have probabilities smaller than this, at 
least in one of the 4 observational constraints. 

5. The evolutionary synthesis models 

In order to derive integrated stellar abundances, we calculate the in- 
tegrated mass of all stars created in a given time interval, and the mean 
abundance reached at that epoch by the gas out of which they form. We 
consider the stellar populations residing at every galactocentric region 
as the superposition of a set of single stellar populations or generations 
each one defined by its age and its metallicity. Spectral index features 
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are calculated with the same method described in Molla Sz Garcia- Vargas 
(2000) by using the Padova group isochroncs. The fitting functions from 
Worthcy (1994) and from Idiart ct al (1997), for assigning the index Fe52 
and Mg2, respectively, to each star are used. We obtain this way the 
radial distributions of these spectral indices for the 19 models of Table 1. 

As before, we compute the by comparing these two radial distri- 
butions with the observations (Beauchamp & Hardy 1997; Molla et al. 
1999). The resulting probabilities are given in Columns (11) and (12) of 
Table 1. Only 6 of the 19 models, marked with an * in this same table, 
fit these new constraints with probabilities larger than 97.5 %. 

6. Conclusions 

We have computed 500 models for the galaxy NGC 4303 with dif- 
ferent input parameters. By using the usual goodness-of-fit chi-square 
parameter we select those reproducing the present-day data within a 
confidence level of 90%, and we constraint the possible number to 19. 
However, these models do not reproduce equally well the spectral index 
radial distributions. Out of these 19 models, only 6 are also able to 
reproduce the radial distributions of spectral indices. 

It seems that the uniqueness problem associated to chemical evolution 
models is not strong: possible models reduce to 5% of the initial ones 
using only the present day data as constraints. When wc also use the 
spectrophotometric indices, we limit even more the possible evolutionary 
histories, by reducing to a 1% the possible models out of the initial 500. 

We conclude that this technique combining chemical evolution with 
evolutionary synthesis is a very powerful tool to discriminate among 
evolutionary scenarios in galaxies with continuous star formation. We 
propose more observational campaigns in order to obtain the present as 
well as other spectral indices in spiral and irregular galaxies. 
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